ORTHOGONAL POLYNOMIALS AND EXPANSIONS FOR A 
FAMILY OF WEIGHT FUNCTIONS IN TWO VARIABLES 



n YUAN XU 

Abstract . Orthogonal polynomials for a family of weight functions on [— 1 , 1] 2 , 
W„ Al (*,») = [x + J/| 2a+1 |x - y| 2/3+1 (l - x 2 y(l - y 2 )\ 

are studied and shown to be related to the Koornwindcr polynomials defined 
on the region bounded by two lines and a parabola. In the case of 7 = ±1/2, an 
^— I explicit basis of orthogonal polynomials is given in terms of Jacobi polynomials 

I and a closed formula for the reproducing kernel is obtained. The latter is used 

to study the convergence of orthogonal expansions for these weight functions. 

1. Introduction 

i-C Orthogonal polynomials of two variables with respect to a nonnegative weight 

function that has all moments hnite are known to exist ([!]). A basis of orthogonal 
polynomials can be written down, say, in terms of moments, but such a basis is 

i_i often hard to work with. For studying orthogonal polynomials and orthogonal 

expansions, additional structures are often called for. In the case of classical weight 
functions in two variables, for example, an orthogonal basis can be expressed in 

qq terms of classical orthogonal polynomials of one variable. There are, however, not 

\^ many such examples; each additional one is valuable in its own right. 

CN| The purpose of the present paper is to study orthogonal polynomials and orthog- 

^ onal expansions with respect to a family of weight functions defined on [— 1,1] , 

fNJ which includes as a special case 

O (1.1) W a ^(x,y):=\x-y\ 2a + 1 \x + y\ 2 ' 3 + 1 (l-x 2 r(l-y 2 r, 

where a, /?,7 > —1 and a + 7 + 3/2 > and (5 + 7 + 3/2 > 0. In the case 
a = (3 = —1/2, W a ,^,7 is the product Gegenbauer weight functions, for which an 
orthogonal basis is given by product Gagenbauer polynomials. We shall show that 
rS orthogonal polynomials for this family of weight functions can be expressed in terms 

of orthogonal polynomials in one variable when 7 = ±1. Our study starts from a 
realization that it is possible to express the orthogonal polynomials with respect to 
Wo, ,3.7 in terms of the Koornwinder polynomials that are orthogonal with respect 
to the weight function 

(1.2) W a> p n [u,v) := (l-u + v) a {l + u + vf(u 2 - 4w) 7 
defined on the domain Q bounded by two lines and a parabola, 

(1.3) Q := {(u, v) : 1 + u + v > 0, 1 - u + v > 0, u 2 > Av}. 
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Orthogonal polynomials with respect to W a ,p,^ were first studied by Koornwinder 
in [7] , where an orthogonal basis is uniquely defined and shown to consist of eigen- 
functions of two differential operators of order 2 and order 4, respectively. In the 
case of 7 = the orthogonal polynomials can be given in terms of the Jacobi 
polynomials of one variable. Further studies were carried out in [9l[T3]; in particular, 
explicit formula for the orthogonal polynomials were derived and various recursive 
relations were established. The connection to orthogonal polynomials with respect 
to W Qjy 3 j7 is somewhat surprising, but simple in retrospect as can be seen by the 
relation 

W a)/3 , 7 0, y) = 4-iW a ,f,„(2xy, x 2 + y 2 - l)\x 2 - y 2 \. 

As a result of the connection, an orthogonal basis for W Q p ±i can be given explicitly 
in terms of the Jacobi polynomials. 

An explicit orthogonal basis makes it possible to study orthogonal expansions, 
for which however it is essential to have access to the reproducing kernel of the space 
of polynomials of degree at most n in L 2 (W). It turns out that closed forms of 
the reproducing kernels for W a p±i and for W a t p t ±± , respectively, can be given in 
terms of the reproducing kernels of the Jacobi polynomials. This allows us to prove 
several results on the convergence of the orthogonal expansions for these weight 
functions. The results include L p convergence of the partial sum operators and 
sharp estimate of the Lebesgue constants. It is interesting to note that analogous 
result on the L p convergence has not been proven for the product Jacobi series 
on the square (the partial sum is defined in terms of polynomial subspace of total 
order, see Remark 2.1). In fact, as far as we know, W a ^_i appears to be the first 
family of weight functions on the unit square for which a comprehensive study of 
orthogonal expansions is possible. 

The Koornwinder polynomials are derived from the symmetric orthogonal poly- 
nomials with respect to the weight function 

W a ,^(x + y,xy)\x - y\ = (1 - x) a (l + xf{l - y) a (l + yf\x - y\^ +l , 

which are the generalized Jacboi polynomials of BC2 type. The latter are the 
first case of the generalized Jacobi polynomials of BC' n type studied by several 
authors (see, e.g. [2J[H]) and they motivated the Jacobi polynomials of Heckman 
and Opdam [6] associated with root systems. The connection between these BC2 
polynomials and orthogonal polynomials for W Qi( g, 7 appears to be new. 

The paper is organized as follows. The following section is a preliminary, where 
the basic results on orthogonal polynomials are introduced. In Section 3 we recollect 
properties of orthogonal polynomials for W a ,p,-y and establish a closed form formula 
for the reproducing kernel. The orthogonal polynomials for W a p ^ are studied in 
Section 4. The orthogonal expansions are investigated in Section 5. 

2. Preliminary on orthogonal polynomials 

In this short section we recall basics on orthogonal polynomials of one variable 
and two variables, respectively, in two separate subsections. 

2.1. Orthogonal polynomials of one variable. Let w be a nonnegative weight 
function on [—1, 1] that has finite moment of all orders. Throughout this paper we 
denote by p n the orthonormal polynomials of degree n with respect to the weight 
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function w, which are uniquely determined by 

J p n (x)p m (x)w(x)dx = 5 n , m , n, m > 

and 7„ > 0, where 7„ denotes the leading coefficient of the orthogonal polynomial 
p n , that is, p n (x) = j n x n + . . .. Let U n denote the space of polynomials of degree 
at most n in one variable. The reproducing kernel of fc„(u>; ■, •) of II„ is defined by 
the relation 

J P{y)k n (w;x,y)w(y)dy = p(x), Vp e n„. 
The well known Christoffcl-Darboux formula shows that 

/O-n u f \ V"" I \ ( \ 7n Pn+l(x)pn(y) - p n +l(y)Pn(x) 

(2.1) k n (w;x,y) = > Pk{x)p k {y) = . 

fc = x - y 

The Fourier orthogonal expansion of / G L 2 (w) is defined by 



f = ^2 f nPn and f n : ^ / 



f(y)Pn(y)w(y)dy, 



n=0 

where the equality holds in the L 2 (w) sense by the standard Hilbert space theory 
and the fact that polynomials are dense in L 2 (w). The partial sum operator s n f 
of this expansion is given by 

(2.2) s n {w;f,x):=^2f k p k = f(y)k n (w;x,y)w(y)dy, 

where the second equal sign follows from the definition of k n (w;-, •). 
The Jacobi weight function w — w a ^ (w e J) is defined by 

w a ,p(x) := (l-x) a (l + x)P, a,/3>-l. 

The Jacobi polynomials are orthogonal with respect to w ai p and they are given 
explicitly as an hypergeometric function 

P W)( X y)- {a + 1)n ,F. (-n,n + a + p + l_l-x' 
n [X ' y) ~ n \ 2ri \ a + l 

These polynomials satisfy the orthogonal conditions 

c a ,p J Pt P \x)pW\x)w ai p(x)dx = hi a ^S nim , 



where 



T(a + p + 2) , p) = {a + l) n (j3 + a) n (a + P + n + l) 

2»+P+^T(a + l)T(|3 + ^) , n ' n!(a + /3 + 2)„(a + /3 + 2n + l) ' 



We denote the orthonormal Jacobi polynomials by Pn ■ It follows readily that 
p^n'^ix) = (h^'^)~ipi a ' P \x). Furthermore, for w = w aj p, we also write the 
reproducing kernel as fci"'' 3 ^-, •) and the partial sum operator as s£ l Q '' 3 ' (/). 

More generally, a function w is called a generalized Jacobi weight function (w € 
GJ) if it is of the form 

r 

(2.3) w(x) = iP(x)(l-xy°{l + x)~<'-+ 1 Y[\x-Xi[", 7i>-l, 
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if \x\ < 1 and w(x) — if \x\ > 1, where —1 < x\ < . . . < x r < 1 and ip is a positive 
continuous function in [—1, 1] and the modulus of continuity uj of ip satisfies 

—^-dt < oo, 

which holds, in particular, if V> is continuously differentiable. For a class GJ, the 
points afi, . . . , x r are fixed whereas 7i, . . . , 7r are parameters. In the case of 7; = 0, 
1 < i < r and VK 21 ) = 1> u> is an ordinary Jacobi weight function. Orthogonal 
polynomials with respect to w £ GJ are called generalized Jacobi polynomials. 
They share many properties of Jacobi polynomials (see, e.x., [J [TO]), even though 
they do not have explicit formulas in terms of hypergeometric functions. 

2.2. Orthogonal polynomials of two variables. Let W be a nonnegative weight 
function defined on a bounded domain Q C M 2 . We define an inner product 



(2.4) (fi9) w := / f{xi,x 2 )g(xi,X2)W(xi,x 2 )dxidx2 

on the space of polynomials. Let II 2 denote the space of polynomials of (total) 
degree at most n in two variables. A polynomial P E II 2 is called orthogonal if 
(P,Q) = for all Q £ n^_ x . Let V n {W) denote the space of such orthogonal 
polynomials of degree n. Then 



dim V n (W) =ra + l, dimlT 2 



n + 2 
2 



The space V n (W) can have many different bases. We usually index the elements 
of a basis by {Pk, n ■ < k < n}. A basis of V n (W) is called mutually orthogonal if 

(Pk,n,Pj,n)w = h k5k, 3 , < k,j < n, 

and it is called orthonormal if hu = 1, < A; < n. The reproducing kernel 
K n (W; -, •) of n 2 in L 2 (W) is defined uniquely by 

K n (W; x, y)f(y)W(y)dy = f(x), V/ e Tl 2 n , 

where x — (21,2:2) and y = (1/1,2/2)- Let {Pk. m ■ < k < m} be a sequence of 
orthonormal polynomials with respect to W. Then the kernel K n (W; •, •) satisfies 

n m 

(2.5) K n (W;x,y) = J2J2 P ^m( x ) P k,m(v)- 

Since W is bounded, polynomials are dense in L 2 (W). For / 6 L 2 (W), the orthog- 
onal expansion of / is defined by 

00 n „ 

f = X] X] fk, n Pk,n, where f k>n = J f(y)Pk,n(v) W (v) d V- 
The n-th partial sum operator of the above expansion is give by 

n k „ 

(2.6) S n (W;f):=J2Y.h> kP ^= / f(y)K n (W;-,y)W(y)dy 7 
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where the second equal sign follows from ( |2.5[ ). There is an analogue of the 
Christoffel-Darboux formula for this kernel ([Up. 109]) but it still involves a sum- 
mation and is not as useful. For studying convergence of the orthogonal expansions 
beyond L 2 , it is often necessary to have a compact formula for the kernel. 



Remark 2.1. The partial sum S n (W; f) in (2.6) is defined in terms of the polynomial 
space 11^ in total degree. Since p n {x)p m {x) has degree n + m, the partial sum for 
the product weight function, say W at p(x,y) — w a ^(x)w a ^{y), does not have a 
product structure. In fact, there is no compact formula for the kernel K n (x, y) for 
W at p if (a,j3) 7^ (— i, — |). As a consequence, there is little progress on the study 
of orthogonal expansions on the square. 

3. KOORNWINDER ORTHOGONAL POLYNOMIALS 

The definition and the properties of the Koornwinder orthogonal polynomials 
are discussed in the first subsection. A new compact formula for the reproducing 
kernels is given in the second subsection. 

3.1. Orthogonal polynomials. Let w be a nonnegative weight function defined 
on [—1, 1]. For 7 > — 1 define 

(3.1) B 7 (x,y) :=alw(x)w(y)\x-y\^ +1 , (x, y) e [-1, l] 2 , 

where a 7 is a normalization constant such that Jj_ 1 ^ 2 B 7 (x, y)dxdy = 1. 

Since i? 7 is evidently symmetric in x, y, we only need to consider its restriction 
on the triangular domain A defined by A := {(x,y) : — 1 < x < y < 1}. Let SI be 
the image of A under the mapping (x, y) 1— > (it, v ) defined by 

(3.2) u — x + y, v = xy. 

It is easy to see that this mapping is a bijection between A and SI. The domain SI 
is given by 

SI := {(u, v) : 1 + u + v > 0, 1 - u + v > 0, u 2 > Av} 



and it is depicted in Figure 1. This is exactly the domain defined in (1.3) 



Figure 1. Domain SI 



We consider a family of weight functions defined on the domain SI by 
(3.3) W 7 (u, v) = 2alw(x)w(y)(u 2 - 4w) 7 , (u, v) G SI, 



where the variables (x, y) and (u, v) are related by ( 3.2 ) . The Jacobian of the change 
of variables (3.2 1 is given by dudv = \x — y\dxdy. Moreover, u 2 — Av = (x — y) 2 . It 
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follows that 

(3.4) / f(u, v)Wj(u, v)dudv = 2 / f(x + y,xy)Bj(x,y)dxdy 



A 

f(x + y, xy)B 7 (x, y)dxdy, 

[-i,i] 2 

where the second equal sign follows since the integrant is a symmetric function of 
x and y, and [— 1,1] is the union of A and its image under {x,y) H ► (y,x). In 
particular, setting f(x, y) — 1 shows that W 7 is a normalized weight function. 

In the case of w(x) = w a ,p(x), the weight function W 1 becomes W a> p^ in (1.2 1, 
which we restate below, 

(3.5) W a , M {u,v) :=2a a ^^(l-u + v) a (l + u + v) l3 (u 2 - 4:v)~<, (ti.ujeO, 

where the constant a Q ,,3, 7 is given by |13l Lemma 6.1], 

, s r( Q + /3 + 7 +|)r(Q + /3 + 2 7 + 3) 

1 ; aQ '^' 7 ' 22^+2^+47+4 r(a + 1)r(/8 + 1)r(7 + 1)r(a + 7 + | )r(/3 + 7 + |) ' 

which is the two- variable case of the Selberg integral (e.x., (1.1) of [S]). This weight 
function is integrable on Q, if a, /3, 7 > —1 and a + 7 + 3/2 > and (5 + J + 3/2 > 0, 
and we assume that a, /3, 7 satisfy these inequalities from now on. Other examples 
of W 7 include 

e' au (u 2 ~4vy and e^" 2 " 2 "^ 2 - 4w) 7 , a > 0. 

which correspond to the choices of w(x) = e~ ax and w(x) — e~~ ax ~ . 

Let A/" = {(k, n) : < k < n}. In A/" define (j,m) -< (fc,n) if to < n or j < fc 
when m — n. Then the orthogonal polynomials Pj^ that satisfy 

(3.7) P^(u,v) = u n - k v k + a ~ - m ~ j " 4 

(j,m)-<(fe,7i) 

and the orthogonality condition 

,(7) 



(3.8) / P^(w,u)w m ^^W 7 (M,i;)dMdi; = 0, V(j,m) -« (fc,n), 

are uniquely determined, as can be seen by the Gram-Schmidt process. The poly- 
nomials Pf^l are mutually orthogonal. 

When W 1 — Wq^ i7 , we denote these orthogonal polynomials by P£ ' 7 . 

In the case of 7 = ±i, these orthogonal polynomials can be given explicitly, as 
can be easily verified upon using (3.4). Let p n denote the orthogonal polynomial 
of degree n with respect to w. Then an orthonormal basis with respect to W_i is 
given by 

( 3 g x P ( ~i\ u v) = — f Pn W Pk ( yS> +Pn(v)Pk(x), < k < n, 

^\V2p n (x)p n (y), k = n, 

and an orthonormal basis with respect to Wi is given by 



(3.10) P^(u,v) = g ^'W^-^ifeW , < k < 

V 2aJ 2 x-y 
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both families are defined under the mapping (3.2 1. It should be noted that the 
polynomials in (3.9 1 and (3.10) are normalized by their orthonormality, instead of 
by the leading coefficient as in (3.8 1, but they have the same structure as those in 
(3.8), so that the difference is just a constant multiple. In the case of w = W ai /3, 

■ . i . a g — i a 8 - 

the orthogonal polynomials in (3.9) and (3.10) are denoted by P k n ' 2 and P k ' n ' 2 , 

respectively, and they are expressed in terms of Jacobi polynomials p n a '^ ■ 

For W a p -v these orthogonal polynomials were first studied by Koornwinder in 
[7j, see also [5]. The above statements for more general weight function W 1 are 
straightforward extensions and used in |12) for studying Gaussian cubature rules. 
Much more can be said about the orthogonal polynomials P k 'n- They are, for 
example, eigenfunctions of two differential operators of order 2 and order 4, respec- 
tively [7]. Another pair of differential operators were constructed in [13], 

d 2 _ 

dv 2 



(3-11) 



E 



a.fi 



d 2 n , 



1] 



dudv 



d 2 d d 

m— + (/? - a)— + (a + j3 + 2), 



dv 



di 



and they act as raising and lowering operators on the orthogonal polynomials, 



E-'^Pkt 1 ^, v) = (n- k)(n + k + a + f3+ l)P^l\ +L (u, v) 



.q,/3,7+1/ 



E 



+ 



Pkn-i (u,v) = (n-k + 2 1 + l)(n + k + a + l3 + 2 1 + 2)P«>™ (u, v) 



for < k < n — 1 and E"' 13 P"'^ n (u,v) — 0. Together these two operators can 
be used to give a Rodrigues type formula for P k '^' J ([131 (5-1)]) an d they can also 
be used to calculate the L 2 -norms of P k ^ and the coefficients in the recurrence 
relation. 

The polynomials P^'I!' 1 { u i v ) a ^ so satisfy a quadratic transformation formula [131 
Theorem 10.1] given by, for < k < n, 



(3.12) 



n+k p7» 2 > 
r k,n 

-n+k pf- h 
r k,n 



'(2v,u 2 -2v-l), 
\2v,u 2 -2v-l). 



In particular, setting 7 = ±^ and a 4 7 and let s = 2xy and t 
follows that p2 p. 518], 

< 1) = c{pi'; k (x)p^ k (y)+PlU^pl'Uy^- 

pi ' 2 ' 

(3.13) 



1, it 



1 

p 2' 

r k,n 


i7 ( S ,o 


1 

pi ' 


i7 ( S ,o 


1 

P 2 







r 1 [pllk 

\~ X IV,7>7 



„7,7 
.,7,7 



7,7 
n 7.7 



-7.7 
J n+k+ 



l(*)] . 



J 1 rp^ a (« fc (y)-p^(« fc+ 2W] 



where c is a constant proportional to 2 n ~ k . In other words, a basis of orthogonal 
polynomials for ±^7 can be explicitly given in terms of Jacobi polynomials. 

including explicit series expansions and recursive 



For further results on P, 



a,/3,7 



relations, see [7[ l9"l H3]. 

It is worth to mention that the relation (3.4) shows that orthogonal polynomials 
Pfc.„ for W 7 are closely related to the orthogonal polynomials with respect to B 1 
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on [— 1, l] 2 , as seen by 

W 7 (x + y,xy) = B 1 {x,y)\x - y\ 



and (3.4 1. Indeed, if Rk,n is polynomial orthogonal to x y m , if max{i, m} < n, with 
respect to B 1 on [— 1,1] 2 , then Rk, n (x,y) + Rk,n(y,x) is a symmetric polynomial 
orthogonal to x l y m + x m y l , if max{Z, m} < n, for B 1 since -B 7 is symmetric. Hence, 
under the bijection (3.2 1, the polynomial 

Pk,n(u, v) := R k ,n{x, y) + Rk,n(y, x) 



is an orthogonal polynomial with respect to W 1 on J7. Since the mapping (3.2 1 is 
not linear, one needs to be careful about the degree of Pk.n- In the case of w = io a ,/3j 
the symmetric orthogonal polynomials for B 1 are the BC2 type polynomials, the 
precursor of the generalized Jacobi polynomials of BC n type. 

3.2. Reproducing kernel. Recall that K n (Wj; -, •) denotes the reproducing ker- 
nel of II 2 in L 2 (W 7 ), which we shall denote by (•, •) below. In contrast to (2.5 ), 
we derive a closed formula for Kn(-, •) in the case of 7 = ±| in this subsection. 



Theorem 3.1. Let k n (-, •) := k n (w; -, •) be the kernel defined in (2.1). Set 

u:=(ui,u 2 ) = (xx+X2,XiX 2 ) and v := {v x ,v 2 ) = (2/1 + 2/2,2/12/2)- 
Then the reproducing kernel K n 2 (■, •) for W_i is given by 

(3.14) K n 2 (u,v) = -{k n (x l7 y 1 )k n (x2,y2) + k n (x2,yi)k n (x 1 ,y2)}, 

(.-) 

and the reproducing kernel K n 2 (•, •) for Wi is given by 

,,,,, 7^(5)/ x _ 1 k n+1 (x 1 ,y 1 )k n+1 (x 2 ,y2) - k n+1 (x2,yi)k n+1 (x 1 ,y 2 ) 
2 (x 1 -x 2 ){yi -2/2) 



Proof. Denote the right hand side of (3.14| by k n ((xi,x 2 ), (2/1,2/2))- By the defini- 
tion of fc„(-, •) in (2.1 ), for a fixed 2/1, 2/2, we have 

^ n n 

k n ((xi,x 2 ), (2/1,2/2)) = 2^2^2pk(yi)pj(yi) \pk(xi)pj(x 2 ) +Pk{x 2 )Pj(x 1 )] , 

" k=0 3=0 

which shows, upon setting u\ — x% + x 2 , u 2 — X\x 2 , that fe n ((a;i, x 2 ), (1/1, 2/2)) is 
a polynomial of degree n in (^1,2/2)- Hence, if we define if„((ui,U2), (^1,^2)) = 
k n ((xi,x 2 ), (2/1,2/2)) under the mapping (1*1,^2) n- (#1 + £E2,a;ia;2) and (vi,v 2 ) n- 
(2/1 + 2/2, 2/12/2), then is a polynomial of degree n in (mi, U2) and, by symmetry, in 
( w i, ^2)- Thus, we only have to verify the reproducing property. For < j < m < n, 
the reproducing property of k n {-, •) implies immediately 

K n ((ui,u 2 ), (v 1 ,v 2 ))P j m 2 (v 1 ,v 2 )W__i(v 1 ,v 2 )dv 1 dv2 

k n {(xi,X2),(yi,y2))P jm 2 (2/1 + y2,yiy2)w(y 1 )w(y 2 )dy 1 dy 2 

-l,i] 2 

-y| \Pm(xi)Pj(x 2 ) + Pj(Xi)p m (x 2 )] = P j<m 2 (U1,U 2 ), 
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which shows that K n ((u\, U2), (fi, U2)) is the reproducing kernel of II 2 , so that 



(3.14) holds. 



Denote now the right hand side of (3.15 1 by k n ((xi, X2), (j/i, 2/2))- Then, for fixed 
(yii 2/2), it is easy to see that 

t n w \ \ y> y> Pj (Vi )Pk (2/2 ) Pj (xi )pk (x 2 ) - Pj (x 2 )pk (x 1 ) 

Kn((xi,x 2 ), [yi,y2)) — > > , 

— — Vi — V2 a;i — £2 

which shows, since the terms for k = j are zero, that k n ((xi,x 2 ), (2/1, 2/2)) is a poly- 
nomial of degree n in (ui, U2), where Ui = x% + 2^2 and «2 = ^1X2- By symmetry, 
the same holds for (2/1,2/2) with (£1,2:2) fixed. Thus, it remains to prove the repro- 
ducing property, which works similarly as in the case of 7 = —1/2 upon using the 
fact that (j/i — J/2) 2 in Wi{ui,u 2 ) — {yi — y 2 ) 2 w(yi)w(y 2 ) cancels the denominators 

in both k n (', ■) and P/ m (yi + 2/2,2/12/2)- □ 

The closed formula for the reproducing kernel allows us to study the convergence 
of the Fourier orthogonal expansions, which will be discussed in Section 5. 

4. Orthogonal polynomials for weight functions on [— 1, l] 2 
In this section we study orthogonal polynomials for the family of weight func- 



tions defined in (3.3) on [— 1, l] 2 , which are closely related to W 7 . The orthogonal 
polynomials are given in the first subsection, their further properties are in the 
second subsection, and a compact formula for the reproducing kernel is in the third 
subsection. 

4.1. Orthogonal polynomials. Let w be the weight function on [—1,1] and let 
W 7 be the weight function defined in (3.3). We then define 



(4.1) W 1 (x,y):=W J (2xy,x 2 +y'-l)\x 2 -y 2 \ 1 (x,y) € [-1,1] 



which is the quadratic transform that has appeared in (3.12). Let f2 be the domain 



( fL3j ) of W 7 . The fact that (x, y) G [-1, l] 2 implies that (2xy, x 2 + y 2 - 1) G Q is 
shown in the lemma below. In the case that W 1 = W QiJ g i7 , the weight function W 7 
becomes, up to a constant, W a ,p,-y defined in (1.1), which we restate as, 



(4.2) W a ,p, 7 (x, y) := 2a a ,^\x - y\ 2a+1 \x + y\ 2fi+1 {l - x 2 y< {I - y 2 y 

where a, (3, 7 > —1, a + 7 + ^ > — 1 and j3 + 7 + \ > — 1. These conditions on the 
parameters, ensuring the integrability of W Q!( 3 j7 , are the same as those for W a p y 
The weight function VV 7 , and W Qi) g j7 , is normalized as shown below. We define a 
region ST := {(x, y) : — 1 < x < y < 1}. 

Lemma 4.1. The mapping (x, y) <— > (2xy, x 2 + y 2 — 1) is a bijection from Cl* onto 
fi. Furthermore, 

(4.3) / f{u,v)W 1 {u,v)dudv= \ f(2xy,x 2 +y 2 -l)W- ( (x,y)dxdy. 
Jn J [— i,i] 2 

Proof. For (x,y) € [— 1,1] , let us write x = cos 9 and y = cos^, < 0,(f) < tt. 
Then it is easy to see that 

(4.4) 2xy = cos(d - (f>) + cos(6» + 0), x 2 + y 2 - 1 = cos(6> - <f>) cos(6> + </>), 
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from which it follows readily that (2xy, x 2 + y 2 — 1) € Q. For the change of variable 
u — 2xy and v — x 2 + y 2 — 1, we have dudv — 4\x 2 — y 2 \dxdy, from which the stated 
formula follows. □ 

As a more general example, we consider w being a generalized Jacobi weight. 

Lemma 4.2. Let w 6 GJ be defined as in (|2.3|). Then 

I x f u k — yx — ik) — yy — ife 

fc=i 



(4.5) W 7 (x, y) = y)|x - y\ 2 ^ +1 \[\l + t 2 -(x~ t k ) 2 - (y - t k f 



x \x + y\ 2 ^ +1 (l-x 2 r(l-y 2 r, 
where ^ (cos 8 , cos (f>) = ip(cos(0 — 4>))ifj (cos (9 + </))). 

We consider orthogonal polynomials with respect to the inner product 

(4.6) (f,9)w '■= / f(x,y)g(x,y)W 1 (x,y)dxdy. 

J[-X,l] 2 

Let V n (W 7 ) denote the space of orthogonal polynomials of degree n with respect to 
the inner product (■, -) w . It turns out that a basis of V„(W 7 ) can be expressed in 
terms of orthogonal polynomials with respect to W 7 and three other related weight 
functions. Recall that (-,-)w ^ s defined in (2.4|. We define three other weight 
functions 

W^\u, v) := 2t£ cliI) (1 - u + «)(1 + « + u)W 7 (u, «), 

(4.7) W< 1 ' )(«,i;):=2o^ 1 , j(l-u + t;)W y (tt,i;), 
Wf'^u, «) := 2a^ (0tl) (1 + u + v)W 1 {u, v), 

where we define 

(4.8) w {hj) (x) := {l-xy(l + x) j w(x), i,j = 0,1. 

Under the change of variables u — x + y and v — xy, l — u + v = (1 — x)(l — y) and 
1 + u + v = (1 + x)(l + y). The three weight functions in (4.7| are normalized so 
that f^W^'i (u,v)dudv = 1. Clearly these three weight functions are of the same 
type as W 7 . 

In the following theorem, we denote by {P^ ■ < fc < n} an orthonormal 
basis of V n (W 7 ) under (-,-} w . For < fc < n, we further denote by P^' 1,1 , 
-fiTi' 10 '' ^"fcTi' 01 tbe orthonormal polynomials of degree n with respect to (f,g) w 
for = W4 M) , W 7 1,0) , W4 ' 1} , respectively. 

Theorem 4.3. For n = 0, 1, . . ., an orthonormal basis of V2 n (W 7 ) is given by 
lQ ( k,l& V) :=^2(2^ + y 2 - 1), < fc < n, 
»<#i.O»M0 :=^ M V -y 2 )Pii 1 i 1 (2^,x 2 + y 2 - 1), < fc < n- 1, 

and an orthonormal basis o/V2n+i(W 7 ) is given by 

(4 1Q) iQ^L+iO^) == b^\x + y)P$°>\2xy,x 2 +y 2 1), < fc < n, 
aO^n+i^' W) : = 6 7 1,0) (^ - y)PS 4 '°(2.x y ,^ + y 2 - 1), < fc < n, 



where by 
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Proof. These polynomials evidently form a basis if they are orthogonal. By Lemma 
|4.1[ for < j < m and < k < n, 

n(i) \ — /p^ p^\ —A a 



W 7 



and, setting F = PQ-iHt 



(lQMn.aO^L = / (x 2 ~y 2 )F(2xy,x 2 +y 2 ~l)W 1 (x,y)dxdy. 
X 'W, J[-l,l] 2 

The right hand side of the above equation changes sign under the change of variables 

(x,y) H> (y,x), which shows that ( lQ^lw 2 Qj%2rn) = 0- Moreover, since (x 2 — 



y 2 ) 2 W~ f (x,y)dxdy is equal to a constant multiple of W^ 1 ' 1 \u,v)dudv, we see that 



1 nW rSrfi \ _/p(7),i,i P (7).i.i\ r x 

v 2 ( *fe 1 2n' 2t «i,2m/ vv — N^fc.n-l ' r j,m-l / w (i,i) ~ u k,j"n,m- 

Furthermore, setting F = ffcn^fcn' 0,1 > we obtain 

(iQkLiQ^L+i)^, =[ {x + y)G{2xy,x 2 +y 2 -\)W 1 (x,y)dxdy, 

\ ' ' W\ / r_i H2 



[-1,1] 

which is equal to zero since the right hand side changes sign under (x,y) <— > 
(—X, —y). The same proof shows also (iQ^\ n , iQ^ m + lj = 0- Together, we 

have proved the orthogonality of iQ^l n and 2Qk'l n - 

Since (a; — y)(x + y) — x 2 — y 2 changes sign under (x,y) i-» (y,x), the same 
consideration shows that (iQ^\ n +i-> 2< 3j 7 2m+i) ~ 0- Finally, we also have 

n (l) n h) \ _/p(7),0,l p(7),0,l\ _ <r r 

l^fc,2n+l' l^j,2m+l/ w — \^k,n ' r j,m /^(O.l) — u k,jVn,m 

I \ _/p(7),l,0 p(7).l,0\ _r r 

which proves the orthogonality of iQ^, 7 2„ +1 and 2Q^\ n+ i- ^ 



For the weight function W Q ,/3,7, we denote a basis of orthonormal polynomials 
by Pfr'n' 7 m the following theorem. 

Theorem 4.4. For n = 0, 1, . . ., an orthonormal basis of V2 n (VV Q . ( 3. 7 ) is given by 
(4.11) 

i^j^y) ^P^ n (2xy,x 2 +y 2 1), < fc < n, 
2 Q^(x, y) ^'^(a 2 - y 2 )P^J +1 "(2xy, x 2 + y 2 - 1), < k < n 1, 
and an orthonormal basis of V2 n +i(W a ,/3,7) *s given by 

2 Q^ 7 + i(^, y) := - Z/)P fc Q : M ' 7 (2^, + j/ 2 - 1), < fc < n, 

where &£'^ 7 := yja a+i ^ +jn /a a ^ n for i, j = 0, 1. 
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4.2. Special cases and properties. In the case of 7 = ±i, we can derive an 
explicit formula for the basis from (3.9) and (3.10), which takes a particularly 
simple form if we change variables to 

(4.13) x = cos9, y = coscf>, 0<9,<fi<ir. 

Corollary 4.5. For 7 = ±\, th e orthonormal basis defined in (4.9) and (4.10) of 
V„(W ± i) satisfies, under (4.13), explicit formulas upon using the relations 

(4.14) P i k ;J ) (2xy,x 2 +y 2 -l) 



V2 



[p n (cos(6» - 0))p fc (cos(6> + <j>)) + p k {cos(9 - 0))p n (cos(0 + <j>))] 



for < k < n, where the k — n term Pn,n^ is multiplied by the constant y/2; 
furthermore, for < k < n, 

(4.15) P$(2xy,x 2 + y 2 -l) 

/ a^, 1/2 p n (cos(9 - 4>))p k (cos(9 + cf>)) - p k (cos(9 - 4>))p n (cos(9 + cf>)) 



2 a 



1/2 



2 sin 9 sin ( 



Proof. This follows immediately from (4.4), (3.9) and (3.10). 



□ 



In particular, the orthonormal basis for the weight function 

W a)/J ,±i (H,s 2 ) = c\x x - x 2 \ 2a+1 \ Xl + x 2 \ 2f!+1 (l - xj) ±l Hl - x\)^ 

on [—1, l] 2 , where c = 2a a p ± 1 4 ± 5 , can be given in terms of the Jacobi polynomials. 

Proposition 4.6. Let a,f3 > —1. An orthonormal basis o/V2n(W Q/ g _i) is given 
by, for < k < n and < k < n — 1, respectively, 



iQk*2n 2 } (cos 6», cos ( 



2<9fc°2n 2) (C0S6»,C0S( 



1 



p^-«(cos(0-0))pi Q ^(cos(0 + 0)) 
+ Pk a ' p) (cos(9 - 0))p£*«(cos(0 + 0))' 



J,(l,l) 

y o 1 



V\-\ (cos(f-«^))Pfe ; (cos(0 + 



. (a+1,/3+1)/ /„ (a+1,/3+1), /„ . 1 \ \ 

+Pk (cos^ - W)p„-i ; (cos(6» + 0)) 

and an orthonormal basis of V2n+i(W Q ^-i) *s given by, for < k < n, 



2Q { k a 2n+l 2) ( COs9 > COS ' 



u(0,l) 

«A4 



(a + y) - 0)K"' P " i; (cos(0 + 0)) 



K/3+1), 



~72 



+p k a ^ +1 \cos(9 - 0))p<^ +1 >(cos(0 + </>)) 



-(x-y) \ P ^ +1 ^(co S (9 - ^))p[ a+hfl) (co S (9 + 0)) 



+ J3 i a+1 ^ ) (cos(0 - ^))p£* +1 '«(cos(0 + 0))] , 
where whenever k — n the polynomial needs to be multiplied by an additional v2. 
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Similarly, an orthonormal basis for V n (W a p i ) can be given explicitly in terms 
of the Jacobi polynomials upon using (4.15). 

By Theorem 4.3 the orthogonal polynomials for W Q ,/3,7 are expressed in terms 
of orthogonal polynomials for W ai /3,j, which in turn are expressed in terms of the 
symmetric orthogonal polynomials with respect to the weight function 

B a ,p„{x,y) = (1 - + xf{l - y) a {l + yf\x - y\ 2 ^ +1 



on [— 1, l] 2 ; see (3.1 1 with w = w a ^ and the remark at the end of Subsection 3.1. 
Both weight functions W Ql/ g, 7 and B a ^^ are defined on [— 1, l] 2 , and they satisfy 

(4-16) W at _ h ^x,y)=VB in>a (x,y). 

Consequently, there is some kind of automorphism among these orthogonal polyno- 
mials. By Theorem 4.4 symmetric orthogonal polynomials with respect to W 7 _ i a 
are given by, for < k < n, 

3 7,-2. Q /-o™„. ™2 i „.2 



P 



,„ (2xy,x 2 + y 2 -l) and (x + y)P?*' {2xy, x l + y 2 - 1) 



of degree 2n and 2n+ 1, respectively. These are, by ( 4.16[ ), symmetric orthogonal 
polynomials for B atCin , from which we can derive orthogonal polynomials for W a , a ,i 
by a change of variables u = x + y and v = xy as shown in the end of Subsection 



3.1. Since x 2 + y 2 



2v, we see that 



P k,n (>, u-2v- 1) and WiP fe '' n 2 ' (2v, u 2 - 2v - 1) 

are orthogonal polynomials with respect to W a)0tl -y of degree k + n. Comparing the 
leading coefficients by (|3.7[), we conclude that 



k,n 



U P, 



7, 5 . a 



\2v,u 2 
O, u 2 



2v - I) 



2v - I) 



on-fc P«. a >T v ) 
Z ^n-k,n+k\ U ' U )i 



n—kpa,a,-f 



.— fe,'< 



k (u,v). 



These are, however, precisely (3.12). These relations translate to orthogonal poly- 
nomials with respect to yV Qj/ 3 i7 on [—1, l] 2 as follows. Let iQ^'f ' 7 denote the orthog- 
onal polynomials given in (4.9) and ( |4.10[ ) but with P^'f' 1 as the monic orthogonal 
polynomial as in (3.7). 

Proposition 4.7. We have the following quadratic transforms, for < k < |_§ J , 

iQr 1/2 ' Q (cos^cos0) = 2LfJ-^g«-Z M _ 2LfJ+2fc ' 



cos cos 2±*V 



sin 9 sin (j)2Ql' n 1 ^ 2 ' a (cos9, cos 



-2L^J- fc sin^sin^ 



2 X 2V [Ii _a J _ fei2n _ 2L ^;_ 



ij+2fe ( C 



9-0 
2 



Proof. From the quadratic transform formulas satisfied by P^'^'^ (111,112), we have 
K-kl + k(* + V, xy) = 2-^ k P^ a (2xy, x 2 + y 2 1), 



(4.17) 



{x + y^K-^n+k+^ + V^y) = 2- n+k P:: r r a (2xy,x 2 +y 2 1). 



If x = cos 9 and y — cos </>, then it is easy to see that 



x + y = 2cos^cos^, 



xy = cos 



cos 



2 8+4 
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Consequently, the equations for iQ/^J/ 2 '" an d iQ'k-im+i* follow immediately from 
pTf] ), ( |47TT] > and ( |4l2| . For iQ^H^ and aQ^+f we use in addition the 
trigonometric identities cos 6* — cos = 2 sin sin ^±2 anc l cos 2 i^S 
sin sin 0. □ 



cos 



2 



In the case ofa = /3 = — 1/2, the weight function becomes 

W_ i ,_ h ^x,y) = (l-x 2 )"<(l-yy, (x,y)e [-1, l] 2 , 

which is the product Gegenbauer weight function. An orthonormal basis of V„ for 
this weight function is usually given by the product Jacobi polynomials 

p k , n & v) = p [ k n \ x )Pn-t(yl 0<k<n. 

In this case, another basis for can be stated as follows: 

Proposition 4.8. For a — (3 — —1/2, the orthonormal basis for V2„(W_i _i 7 ) 
is given by 



(4.18) 



i (p£# i 7 ^ (y) + p<#> (rf^ (*)) , < k 



V2 



and the orthogonal basis for V2 n _|_i(W_ i _ i 7 ) is given by 



(4.19) 



„(7.7) 



(7,7) 



J7.7)/ 



< n, 

< fc < n - 1, 

< k < n, 
< k < n. 



Proof. The orthogonality of these polynomials follows from (3.13) and Theorem 
|4.4| They can also be verified directly. 



□ 



Finally, let us mention that under the change of variables u — 2xy and v = 



y — 1, the operator E_ in (3.11) becomes 



(4.20) 



?<*,0 ._ 



1 d 2 



1 



2 dxdy 2(x 2 — y 2 ) 



((a + /3 + l)x + (a - j3)y) 



Ox 



((a - P)x + (a + /3 + a)y) 



d_ 
9y 



which has a simple form for the second order derivatives, so that, by (4.11), 



£°'"lQZ$?(x, y) = -(n-k)(n + k + a + f3 + l^Qf^-V^ v)- 
The operator does not, however, act on 2Qt : 2n m ^ ne same manner. As the 



operator E^ has the same second order derivatives as that of E"' p , we can also 
have an that has simple second order derivatives and act on \Q°^2n according 
to |3Tll. 



a,fl 
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4.3. Reproducing kernel. We express the reproducing kernel for K n (W 7 ; ■, ■), 



which we denote by /cl 7 ^-, •) below, in terms of the reproducing kernel Kn\-, •) 
defined in Subsection 3.2. For Wy defined in (4.7) with i,j — 0, 1, we denote by 
K^'^i; •) the reproducing kernel K n (wJf' j) ; •, •)■ 

Theorem 4.9. For x = [xi,x 2 ), y = (1/1,2/2)7 define 

(4.21) s = (s 1 ,s 2 ) = (2x 1 x 2 ,x 2 1 +x 2 2 -l), t=(t 1 ,t 2 ) = (2y 1 y 2 ,y 2 + y 2 .-l). 
Then the reproducing kernel JC^\-, •) for W 7 is given by 

(4.22) K$\x, y) = (s, t) + d^ (x 2 - x 2 )(y 2 - y\)K^{s, t) 



(7)/ 



L L*J- 

+ <f (1 , 0) ( Xl + x 2 )( yi + y 2 )K^(s, t) 

+ <Cco,i) (^i ~ x 2)(yi - y2)K ( ^^(s, t), 



whe 



4*1 — n l 



i,j)/o>1 for i,j = 0,1. 



Pro of. W e consider /C^ (x,y). By the definition of K^\-, ■) as in (2.5) and Theo- 
it follows readily that IC^J (x, y) belongs to IIa n as a function of either x 



4.3 



or y. To see that it reproduces polynomials in n^, we verify 

(K%M,iQ%hy»=iQt!L 0<k<m<2n, i = 1,2, 



using ( |4.3[ ) and ( |4.6[ ). For iQ^l m , this follows immediately from ( |4.3[ ) and the repro- 
ducing property of i^« 7 ^ , since among the four terms in the right hand side of (3.14), 
only the first term has a non-zero inner product with iQ^l m by orthogonality. For 
2 Q\ 2m , we use ( 4.10 ) and, in addition, d 7 (1 1} (x\~ x 2 ) 2 W 7 (xi, x 2 ) = Wy'(ui,u 2 ). 
The other two cases, iQk' 2 m+i W1 th i — 1,2, work out similarly. □ 

In the case of W a ./3,-y, the formula for /C"'^ ,7 (-, •) takes the form 
(4.23) JC^(x,y) = Kffi(s,t) + d^xj - x 2 )(y 2 - y 2 2 )K^^ '(s,t) 



+ <4°fc(*i + x 2 )( yi + y 2 )K a l ^( S ,t) 



a, 18+1,7/ 



+ ^)(yi - y 2 )K^{s, t), 



a+1,/9,7 / 



where d 



(id) 



3.1 



i a +i,fi+3ft/ a <x,fi,~i f° r = 0, 1. In the case of 7 = ±J, we can then 
to deduce closed formulas for the reproducing kernel ICn^^ 2 (•, •)> 



use Theorem 
which take simpler forms in the variables 

(xi,x 2 ) = (cos#i,cos#2) and (2/1,2/2) = (cos</>i,cos0 2 )- 
Indeed, using the relation (4.4) and (s,t) in (4.21), it follows from ( 3.14| ) that 

(4.24) lC P, ~ k (s,t) 

= i [fc^^cos^ - 2 ), 008(0! - fo))*^ (cosfa + 2 ), cos(0i + (j> 2 )) 

+A£^(oos(0 x - 2 ), cos(0! + fc,))A£^(cos(0i + 2 ),cos(0i - 0a))] , 
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and, since cos(#i — 9 2 ) — cos(#i + 6> 2 ) = 2sin#i sin# 2 , 

(4.25) Kn P '^{s,t) = -(sm0isin0 2 sin<Miii<fe)~ 1 
8 

fc"+i(cos(0i - 6» 2 ),cos(0i - 2 ))fc"+i(cos(6'i + 6» 2 ),cos(>i + </> 2 )) 

-*£&(cos(0i - 2 ), 008(0! + 0a))A£&(cOS(0i + 2 ), 008(0! - 2 )) 



Substituting (4.24) and (4.251 into (4.231 gives a compact formula of K n m 



terms of the reproducing kernels of Jacobi polynomials. 

In the case of a = ft = —1/2, the weight function W nj3 _i is the product 
Chebyshev weight 

(xi,x 2 )= - , \ , {x 1 ,x 2 ) e[-i,i} 2 - 



w_ 



2 > 2 ' 2 



Even in this case, the formula (4.23) is new. Previously, another closed formula for 



the kernel K n (Wo', •, •) was given in [17) . Our new formula, however, is more easily 
adopted for studying convergence of Fourier orthogonal expansions as seen in our 
next section. 



5. Fourier Orthogonal Expansions 

In this section, we study orthogonal expansions for both W_ 1 on fi and W_ 1 on 
[—1,1]. The results include both LP convergence and the uniform convergence. The 
LP convergence will be established for W_i and W_i associated with the gener- 



alized Jacobi weight defined in (2.3). The uniform convergence will be established 



for W_i and W_i associated with the Jacobi weight. 

22 

We are mainly interested in the case of W_i, which lives on the square [— 1, l] 2 . 
The study of the L p convergence of the Fourier orthogonal expansion on the square 
has been lagging behind, perhaps unexpected, the study on the triangle and on the 
disk. In fact, the LP convergence for the product Jacobi weight on the square has 
not been established. One reason is the lack of a useable formula for the reproducing 
kernel, which, as we explained in Remark 2.1, does not have a product structure. 
Given this background, our result (see Subsection 5.3) is somewhat surprising, as 
it shows that the case of the weight function 



W. 



(x, y) = \x- y\ 2a+1 \x + y\W +1 (l - x 2 )^ (1 - y 2 y 



or more general W_i in (4.5), can be worked out so much easier than that of the 



product Jacobi weight function on the square. 

To get to W_ 1 on [—1,1], we need to deal with W_i on SI first, which in return 
relies on results on w € GJ . In our first subsection we prove some results for the 
generalized Jacobi series of one variable, which are then used to study orthogonal 
expansions for W_ 1 in the second subsection. The results for W_ 1 are presented in 
the third subsection. Throughout this section, the constant c will denote a generic 
constant, its value may change from line to line, and we denote the ordinary Jacobi 
weight function by J a ,p, that is, 



J a .p{x) := (l-x) a {\ + xY 



-1 < x < 1. 
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5.1. Orthogonal expansions in generalized Jacobi polynomials. Let w be a 

generalized Jacobi weight, w € GJ, defined in (2.3 1. For 1 < p < oo, the LP norm 
of / € W(w, [-1, 1]) is defined by 



l/lltu.P 



i \ Vp 

\f{y)\ p w{y)dy 



1 < p < oo, 



and for p — oo, we replace the LP space by C[— 1,1] with the uniform norm 
||/IU,oo = ||/||oo- For 1 <p < oo let q be defined by l/p + l/q = 1. 

Recall the partial sum operator s n (w; f), see (2.2), of the orthogonal expansion. 
For proving the mean convergence of S n (W_ i ; /), we need to study the convergence 
of a family of operators closely related to s n (w), w G GJ. For i, j > 0, we define 

(5.1) s l ^{w]f,x):=J^j_{x)J f(y)k n {J i!j w;x,y)w(y)J^i(y)dy. 

Evidently, s n (w) = s^°(w). We shall show that these operators have the same 
convergence behavior as that of s n (w; /). 

Standard Hilbert space theory shows that s n (w; /) converges to / in L 2 (w) norm. 
The following theorem gives the convergence of s]^ (w; /) in LP space. 

Theorem 5.1. Let w,u,v G GJ. Then for 1 < p < oo, 

(5.2) \\s l n \wJ)u\\ w , p <c\\fv\\ WtPl i,j>0. 
for every f such that \\fv\\ WtP < oo if and only if 

vPw&L 1 , u p (wJi i )"5wj e L 1 , 

(5.3) v^weL 1 , v- q (wJi i)~iw e L\ 

v 2 ' 2 

and u(x) < v(x), x G ( — 1, !)• 



Ln particular, (5.2) implies that \\{s n (w; f — f)u\\ w _ p — > w/ien n — > oo for every f 
such that HZ^H-u^p < oo. 

Proof. For i = j = 0, this result was proved in [16] (for various earlier results, see 
[TT] and the references therein). We show that the general case of s]^{w; /) can be 
deduced from the case i = j = 0. The operators s\^{w) can be expressed in terms 
of the partial sums of Jacobi series. Let us define 

f itj (x) :=/(i)/ J 4 i(x). 



Then directly from its definition (5.1), we see that 

(5.4) sll 3 (w;f,x) = Ji i{x) s n {Jijw; fij,x). 



The inequality (5.2) is easily seen to be equivalent to, using (5.4), 
( 5 - 5 ) \\s n (Ji.jw; f)u i . j \\ Ji jWtP < c\\fvij\\ 

if we define Uij and Vij by 

Ui,j(y) ■■= Jii(y)(Jij(y))~*u(y) and v^(y) := Jii(y)(Ji,j(y))~*v(y). 



The inequality ( |5.5[ ) holds, by the result for i = j = 0, under the condition (5.3) 
with w replaced by Jijw, u and v replaced by Uij and Vij. We now verify that 
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these conditions hold under (5.3 1. The condition Uij(y) < Vij(y) holds evidently 
under u{y) < v{y) of (5.3). A quick computation shows that 



so that, using Wi.j(y) < c, both are L 1 functions under (5.3 1. A similar computation 
shows that 

U p itj {J i+ ^ j+ iw)-ij id w = u v {Jx t xw)-*w, 

Vi](J i+ i J+ iw)-ij itj w = v- q (Jiiw)-^w. 

Since the right hand sides of the these two expressions are exactly those appeared in 
( 5.3 ) , all conditions under which (5.5| hold are verified under ( 5.3 1 . This establishes 
jgg . □ 

The special case u(x) = v(x) = 1 and w is a multiple of the Jacobi weight is 
stated below as a corollary, in which the conditions in (5.3 1 are simplified to (5.6) 
below. 

Corollary 5.2. Let w(x) — tp(x)(l — x) a (l + x)@ , where a, j3 > —1 and %jj is as in 
(2.3), and let i,j > 0. Then 

||^'K/)IU >f ,<c||/|U, p 

for every f £ L p (w, [—1, 1]) if and only if 
( 5 - 6 ) 2-- r 2 m ■ - <p<2 



2max{a,^} + 3 2 max{a, /?} + 1 ' 



The Theorem 5.1 settles the mean convergence of s]^{w] f) for 1 < p < oo. For 
the cases of p = 1 or p = oo, it is easily seen that 

(5.7) Ha^'HIL = l|s„HIU,i 

= max Ji i(x) \K(J%,jW;x,y)\Ji i(y)w(y)dy. 

xS[— 1,1] 2 ' 2 J — \ 2 ' 2 

In fact, the proof in the case that i = j = is classical and it carries over just as 
well in the general case. We shall determine the order of s % ^ for the classical Jacobi 
weight w = J a ,p and denote, for simplicity, 

By definition, s^'^ = s^"'^' 0,0 is the partial sum of the classical Jacobi series. 
The quantity ||sn ||oo> sometimes called the Lebesgue constant, determines the 
convergence behavior of f when p — 1 and p = oo. 

The asymptotic order of ||s^"' /3 '' Hoc is usually determined by using the convolution 
structure of the Jacobi series, which shows that the maximum 

Halloo = max / \k^(x,y)\J a , p {y)dy 

x£[-l,l] J -l 

is attained at the point x — 1. The same scheme, however, does not apply to 
l '^||oo if either i > or j > because of the factor (1 — x)^ 2 {l + x)' J l 2 in 
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front. Nevertheless, the result still holds and it can be proved by using a sharp 
estimate of the kernel function of s^*'^' 1 ' 3 (/) given by 



COSW, COS ( 



= (1 - x)i(l + x)i(l - y)Hl + y)h^ +j \x,y) 
(sin f sin f ) 4 (cos § cos t^k^^ (cos 9, cos <j>). 



Evidently, '°Gj •) = k^'^ (•, •). It turns out that the kernels in this family 

have the same upper estimate. 

Lemma 5.3. Let a, f3 > — 1/2 and i, j > 0. Then 

(5.8) Ifc^^'^cos^cos^l 



< c 



(sin | sin | + n 1 1 



-2\-a-i/„ Q £ __„ 



(cosf cosf +n" 1 |6» 



1^-01 



Proof. In the case of = (0,0), the estimate is derived from [3l Theorem 2.7] 
by setting S = and d = 1, changing variable from [0, 1] x [0, 1] to [—1, 1] x [—1, 1] 
and then to (cos 9, cos (ft) in [0,7r] x [0,7r], and applying elementary trigonometric 
identities. There are two terms in the estimate in Theorem 2.7] but it is not 
hard to see that the second term is dominated by the first one when 5 — and 
d = l. 

and 



For the general case of i,j > 0, applying (5.8| to fci a+4,,3+ ' 7 ' (cos 9 : cos i 



using the factors (sin | sin |) 1 and (cos § cos |) J to reduce the exponent —(a + i + 
1/2) to -(a + 1/2) and the exponent -(/3 + j + 1/2) to -(/3 + 1/2), we see that 

fcn ^ (cos 0, cos <^>) has exactly the same upper bound as the right hand side of 

(5.8). Consequently, ||sn /3 ^' IJ ||oo has the same upper bound. 



□ 



Theorem 5.4. Let a, j3 > -1/2 and i,j > 0. Then 

r n m<«{a,/3}+l/2 ! max {a, /?} > -1/2, 



(5.9) 



„(Qi,j9),i,j| 



0(1 



logn, 



max{a,/3} = -1/2. 



Proof. The case i = j = can be deduced from the convolution structure of the 
Jacobi series and the Lebesgue function at x = 1 ([HI Section 9.41]). Our proof 
below uses the kernel estimate in (5.8 1 and works for the general case of i,j > 0. 



For i, j > 0, we can rewrite the norm in (5.7) as 

-■l 



\st 0)Aj \\oo = max 



6[0,1] J-l 



k^^(x,y) w a ,e(y)dy 

(sin|) 2Q+1 (cosf) 2 ' 3+1 < 



COS 0, COS < 



By symmetry, we consider only < 9 < ir/2. If ~ < (ft < 7r, then \9 ~ cft\ > 7r/4, 
so that, by the estimate of the kernel, ^"'''''(cos^cos^)! < c as a + | > and 
j3 + h > 0, so that the integral over [37r/4,7r] is bounded. On the other hand, if 
< (ft < 37r/4, then cos | cos | ~ 1, so that the estimate of the kernel shows that 

,.(«,«/ (smfsinf+n-^-^l+n- 2 )-"-^ 
(^^(cos^cos^)! < c * — — — . 
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The case a = — 1/2 is easier, let us assume a > —1/2. We then divided the integral 
into three terms, 

^3?r/4 pB/2 r-38/2 ^3tt/4 

■■■d(f>= ■■■d(f>+ ■ ■ ■ d<f> + / ■ • • dxj>. 

JO Je/2 J30/2 

Using the estimate of the kernel, these three integrals can be shown to be bounded 
by cn a+1 l 2 by making use of the following facts: For the integral over [0,0/2], 
\4>-8\ ~ 0; for the integral over [9/2, 30/2], 0-0; for the integral over [30/2, 3tt/4], 
|0 — <t>\ > 0/3- We leave the details to interested readers. □ 

It should be remarked that the asymptotic order of ||sn ^||oo is established for 
a,/3 > —1. The reason that we assume a,f3 > —1/2 lies in the fact that the 
estimate of the kernel in [3] was established under this assumption. We expect that 
the result extends to 

ll^'^'lloo =0(l)logn, max{a,/3} < -1/2. 

In fact, our proof already shows this estimate if both i,j > 1. Only the cases i = 
or j = remain. What is of interest is to extend the kernel estimate (5.8|, or in 
some modified form, to the range max{a, j3} < —1/2. 

The reason that we restrict to the classical Jacobi weight function in our estimate 
of || s^(w)|| 00 is again the lack of a pointwise estimate for the kernel function. 

5.2. Orthogonal expansions for W_l on fl. Let W be a weight function defined 
on C M 2 . We denote by L P (W) the L p space of functions for which the norm 
is finite, where 

\\f\\w, P := (J \f(x)fW(x)dx^j ^ , 1 < p < oo, 

and, for p = oo, we replace L P (W) by C(f2), the space of continuous functions on 
n with the uniform norm ||/||vy,oo = ||/||oo on f2. 

In this subsection we consider the convergence of S n (W_i ; /) with, see (3.3), 

W_i(u,v) = 2a w w(x)w(y)\u 2 - 4u|~2, w € GJ, 

where u = x + u, v = xy and normalization constant. Because of what we 

will need in the following subsection, we also define a family of operators as follows: 
For i, j > 0, 



(5.10) S^(W_x;f,x):= / f(y)K^(W_x;x,y)W_i(y)dy, 

Jn 

where 

(5.11) K^(W_x;x >y ) := J* 3 (x)J* dy)K n (JtjWi;x,y) 

2 2 , 2 2 ' 2 2 

and 

J* ap {x) := (l- Xl + x 2 ) a {l + Xl + x 2 f. 

Evidently, S®'°(W_i ; f) — S n (W_i ; /). We will show that operators in this family 
have the same convergence behavior. 

Let u, v be the generalized Jacobi weight functions. We define weight functions 
U and V by 

(5.12) U(x + y,xy) — u{x)u(y) and V{x + y, xy) = v(x)v(y) 
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for — 1 < x < y < 1. The functions U and V are well defined on f2. 
Theorem 5.5. Let w,u,v £ GJ and i,j > 0. Then for 1 < p < oo, 
(5-13) \\S^(W_i;f)U\\ w _ ltP <c\\fV\\ w _ ltP 



for all f such that \\fV\\w i,p < 00 if u,v,w satisfy the conditions in (5.3). In 

- 2 ' 

particular, (5.13) implies that \\(S n (W_ i ; /) — iOV^Iw x , p — >• w/ien n — > oo /or 

1 r 2 ' 

all f such that ||/V^|w i , p < oo. 

~ 2 ' 

Proof. For / defined on 57 we define F(x±,X2) '■— f( x i + x ii x i x 2) f° r (^1,^2) € 
[— 1,1] 2 . We first consider the case i = j = 0. Let W*(x) := u;(a;i)u>(;r2)- By 
Theorem 13. 11 we have 



<S n (W 1 ;/,xi + x 2 ,xiX2) = - [S n ^(W*]F,xi,x 2 ) + S n . n {W*;x 2 ,xi)\ , 
2 A 

where S n n (W*', F) denotes the partial sum of the product generalized Jacobi series 
that has degree n in each of the variables X\ and x 2 . By definition, 



S n>n (W*;F,x u x 2 ) 



i-iJ-i 

Similarly, define U*(y x ,y 2 ) := u{y{)u{y 2 ) and V*(y 1 ,y 2 ) := v(yi)v(y 2 ). Then, 
applying (3.4) with 7 = — ~, we obtain upon using the definition of W_i, 

\\S n (W_i_; f)U\\ w 1>p < \\S nin (W*;F)U*\\ w *, p . 
2 2 

where the norm of the right hand side is taken over [— 1, l] 2 against the weight 
function W* . Setting t n (w; x±,y 2 ) := s n (w;F(-,y 2 ),Xi), we can write 

S n , n (W*;F, xi,x 2 ) = s n (w;t n (w;xi, -),x 2 ). 

Thus, the product nature of S n)Tl (W*) allows us to apply Theorem 
conclude that 

\\S n , n (W*;F)U*\\ w * iP < c\\FV*\\w*, P = 4fV\\ w _ ltP , 



5.1 



twice to 



where the equality follows from (3.4). This completes the proof when (i, j) = (0,0). 

The above proof carries over to the case i,j > 0. Indeed, since (xi,x 2 ) i-> 
(xi + x 2 ,xix 2 ) sends 1 ± x\ + x 2 to (1 ± xi)(l ± x 2 ), by Theorem 3.1 we have 

S n (W_i;f,x 1 +x 2 ,x 1 x 2 ) = l - [S^ n {W*;F, Xl ,x 2 ) + S^ n {W*;F,x 2 ,x x )} , 



where S% 3 n (W*; F) can be expressed in terms of (w) at (5.4) exactly as S n , n (W*; F) 
is expressed in terms of s n (w). Consequently, the same proof applies and the desired 
result follows from Theorem 15. II □ 



In the case of p = 1 or 00, we restrict to the case of W_i = W a p _i/ 2 . 
Theorem 5.6. Let a, (3 > — \ and i,j > 0. Then 
(5.14) ||^(I4^,-i)IU = WStlKW^-^Ww^,! 

f n max{ a ,/i) + i/2) ] m&x{a,f3} > -1/2, 



0(1) 



log n, 



max{a,/3} = —1/2. 
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Proof. The standard argument shows that 

11^(^,0,-1)11=0 - \\s^(w aJ ^ k )\\ w ^_ iA 

= m&x(l - Xl + x 2 )Hl + x 1 + x 2 )i / \K n {W a+ii p +j _i-,x,y)\W a+ i 0+ i _dy)dy. 



Applying (3.4), ( 3.14 ) and (5.4), it follows readily that 

ll^'(^„i)|U< (|| S ( Q ' 

so that the stated result follows from (5.9). 



a 



The case of S n (Wi) is harder to work with because of the denominator x — y 
in its kernel (3.15). The proof of Theorem 5.5 could carry over only if we modified 
the definitions of U and V to include a power oi \x — y\. Since it adds little to our 
understanding, we shall not pursuit it any further. 

5.3. Orthogonal expansions for W_i on [—1, l] 2 . We now consider the conver- 
gence of the Fourier orthogonal expansions with respect to W_i on [— 1, l] 2 . 

Let u, v £ GJ be the generalized Jacobi weights. We define the weight functions 
U and V on [-1, l] 2 by 

U (cos 9 , cos cj)) :=u(cos(9 — 0))u(cos(6> + 0)), 
V(cos#,cos0) :=i>(cos(6> — <j)))v(cos(9 + <j))). 

Theorem 5.7. Let w,u,v £ GJ. Then for 1 < p < oo, 
(5-15) \\S n (W_i;f)U\\ W l , p < c\\fV\\ W x<p 



for all f such that ||/V||w i ,p < 00 if u,v,w satisfy the conditions in (5.3). In 

~ 2 

particular, (5.15) implies that \\(S n (W_i; /) — /)V||yv i ,p ~ ► when n — > oo for 
all f such that ||/V||w i ,p < 00 ■ 

~ 2 ' 

Proof. We consider the case <S2n(W_±;/). For a given / on [— 1, l] 2 , we define /* 

by 

f*(2x 1 x 2 ,xj +x% - 1) = f[x\,x 2 ). 

By (4.4), /* is well defined on 17. If t\ = 2y\y 2 and t 2 — y\ + y 2 — 1, then it is easy 
to see that 1 — t\ + t 2 = (yi — y 2 ) 2 and 1 + t\ + t 2 = [y\ + y 2 ) 2 . Hence, by (4.22) 
and (5.11 1, we obtain that 

/C n (W_ i ; x, y) =K l n 1 (W_ i ;s,t) + d^K ^ (W_x ; s, t) 

+ d ltl K 1 l '^ i (W_ i -,s,t)+d 1 , 1 K[^ i _ 1 (W_i;s,t). 
Consequently, by ( |4.3| , it follows that 

(5.16) S 2n (W_ i;f,x) = S n (W_ h ;f*,s) + S 1 n l 1 {W_ i ; /*, s) 

+ (W_ , ; /* , a) + S^iW. i;f*,s), 



where s = (2xix 2 , xf + x?,— 1). Recall the definition of U and V in (5.12 1. A simple 
computation via elementary trigonometric identities shows that 

U(xi, x 2 ) — U(2xix 2 , x\ + x\ — 1) and V{x\, x 2 ) = V{2x\X 2 , x\ + x\ — 1). 
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Consequently, applying (4.3 1 again, we see that 

ll«53n(W_i;/)W|| w 1 , P <\\S n (W_ L ,r)U\\ w 1J , + \\SW 1 QV-i;r)U\\ w liP 

2 2 2 2 2 2 

+ \\S n 'l 1 (W_ l ;r)U\\ W rv + WS^WrlDUWw ltP 

2 2 2 2 

<c\\f*V\\w_,, P = c\\fV\\y V _ k , P 



by Theorem 5.5 with i,j = 0, 1. The proof for S2n+i{W_if) follows analogously. 

□ 

Corollary 5.8. Letip be a positive, continuously differentiable on [—1, 1] and define 
*S>(x,y) by ^(cos 9, cos</>) = ^/>(cos(0 — <f>), cos(9 + </>)). Then for 

Wcpfay) := *(x, y)\x - y\ 2a+1 \x + y\^ +1 (l - x 2 )-i(l - y 2 yi, 
where a, j3 > —1. 

||<S , n(Wo, J a;/)||wa, J 5,p ^ c ll/l!w Qi/3 ,p 
for all f such that \\f\\w a $, P < 00 */ (5.6) holds. 

Even for the case a = j3 = —1/2, this corollary is new. In the case of a = (3 = — | 
and i/j(x) = 1, i.e., the product Chebyshev weight function, this was established in 
|18) by identifying the orthogonal expansion with the double Fourier series and S n f 
with the l\ partial sum of the double Fourier series, then applying several results 
in the Fourier analysis. It is worth mentioning that no analogous results are known 
for the product Jacobi weight J a ,p{x)J a ,p{y) on the square. 

For the norm with p = 1 or oo, we go back to the weight function associated 
with the Jacobi weight. 

Theorem 5.9. Let a, j3 > -1/2. Then 

(5.17) llSnfW^.^Hoo = \\Sn{W a ^)\\ W ^_ h + 

max {a,/3} > -1/2, 



0(1) 



log n, max{a, /3} = — 1/2. 



Proof. This is an immediate consequence of (5.16) and Theorem 5.9 □ 
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